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Abstract 



We demonstrate on simple examples that oscillatory behaviour of moments 
of multiplicity distributions P(n) observed in e + e~ annihilations, in hadronic 
pp collisions and in collisions on nuclei, p + A, is to a large extend caused by 
the experimental artifact of measuring only limited range of P(n). In partic- 
ular we show that by applying a suitable universal cutt-off procedure to the 
measured P{n) one gets for reactions mentioned before oscillations of similar 
magnitude. The location of zeros of oscillations as a function of the rank of 
moments and their shapes remain, however, distinctively different for different 
types of reactions considered. This applies to some extend also to collisions 
of nuclei, which otherwise follow their own pattern of behaviour. 

PACS numbers: 13.65,+i 13.65.-7 24.60.-k 
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The problem of the possible physical origin and information content of oscilla- 
tions in the cumulant moments of the corresponding multiplicity distributions P(n) 
started with QCD calculation of the respective generating functional. It turned out 
that the resulting cumulant moments oscillate as a function of their rank in the 
way depending on the QCD parameters used This finding was confirmed by the 
analysis of e + e~ and hadronic pp data, which showed that, indeed, the g-th rank 
normalized cumulant moment of observed negatively charged multiplicity distribu- 
tions oscillates irregularly around zero with increasing q (the minimum points being 
located around q ~ 5). There was therefore a hope that analysis of these oscillations 
can then prove a crucial test for the QCD 0. 

These expectations were, however, soon confronted with observations that the 
same data can be equally well described by a more phenomenological methods H [| 
based on the solutions of stochastic processes || , for example by the negative bino- 
mial distribution (NBD) |6j (in its truncated version) or by the modified negative 
binomial distribution (MNBD) || [F). The interesting finding was that the NBD and 
MNBD differ distinctively in this context in the following sense: for untruncated 
multiplicity distributions the q-th rank normalized cumulant moment of the NBD 
is always positive and decreases monotically with increasing q, whereas for MNBD 
it can oscillate in the way depending on the choice of parameters. Therefore, in 
this approach, the behavior of cumulant moments obtained from experimental data 
seemed to provide a new constraint on models of multiplicity distributions. In par- 
ticular it was shown in H || that cumulant moments of negatively charged particles 
in e + e~ collisions can be described by both the truncated NBD and (truncated or 
not) MNBD (which performs better in the e + e~ case). 

The observation that truncation of the NBD makes the corresponding moments 
oscillate has confirmed the statement made before in Ref. ||. It was said there 
that important (if not exclusive) factor leading to oscillations of moments is the 
experimental fact of necessary truncation of the observed P(n) at some maximal 
multiplicity n max . The observed differences between results from different reactions 
seem to reflect therefore only the level of this truncation. 

Cumulant moments obtained from the hA experimental data also oscillate with 
magnitude which is much bigger than that observed in hh collisions ||. The specific 
feature of these oscillation is that they can be atributed not only to truncation of 
P(n) but also to the fact that the number of elementary collisions in the hA reac- 
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tion is necessary limited by the geometry of collision ||10|| . This is, in fact, a kind of 
truncation as well, but this time it is caused by the geometry of the collision rather 
than by the experimental setup. This new geometrical factor should be therefore 



even more important in heavy ion collisions [11, 12 



The above discussion clearly shows that information content of the oscillation 
phenomenon remains still unclear and subject to debate. The aim of our note is thus 
to shed a new light on this problem by discussing a couple of simple but illustrative 
numerical examples of oscillations in e + e~ annihilations, hadronic pp collisions and 
collisions involving nuclei, pA. We shall not attempt here a fit to experimental data 
because this was done already in the relevant works quoted here. Our intention was 
rather to use the existing experience on this subject (especially that contained in 
H [|, [H|) in order to demonstrate a possible universality existing in the e + e~ , pp 
and pA data on oscillations of moments. We shall also address, albeit only shortly, 
heavy ion collisions, AB, in this context. 

We shall use, as is usually done, the following moments of the multiplicity dis- 
tribution P(n) (cf. 13): 

K 



q F 



where 



with k q and f g being the usual cumulant and factorial moments of rank q of P(n). 
What is observed experimentally is the fact that they oscillate and that these os- 
cillations differ substantially depending on the type of reaction and do it in two 

ways: 

• their amplitudes vary, increasing from the value of 1CT 3 for e + e~ annihilations, 
via the value of 10 -2 for pp collisions up to the value of 10 -1 for pA and heavy 
ion (AB) reactions; 

• their shapes are different with frequency of oscillations in q being highest for 
e + e~ reactions. 

As was said before, the main cause of these oscillation is supposed to be experimental 
truncation of the corresponding P(n). In order to compare results of such truncation 
for different reactions we propose to use a universal variable u defined in the following 

way: 

u = . (3) 

cr-n 
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This variable measures distance of the cut-off point, n max , from the (specific for the 
process under consideration) mean multiplicity (n). It does this in terms of stan- 
dard deviation <j n (which is obtained from the same P(n)). In this way it allows to 
compare results from different reactions by providing a kind of natural and universal 
measure for terminating P(n) under considerations at some value of n max . 

In Fig. 1 we show examples of H q moments calculated for e + e~, pp and pA 
reactions for two different choices of the values of cut-off parameter u. In each case 
we have used identical multiplicity distributions (and all other relevant parameters) 
as those used in Refs. ]3], £|, |10J| when describing the same reactions. They were then 
cut-off for each reaction considered at the same value of the variable u defined in 
eq. (0) and from them the corresponding moments H q were calculated. As can be 
seen, cutting off multiplicity distributions P{n) (calculated for different reactions) 
at the same values of u results in comparable values of amplitudes of observed oscil- 
lations. Although they are still not identical, the previously mentioned differences 
in amplitudes are enormously reduced, being now of the same order of magnitude. 
This feature apparently does not depend on the actual value of variable u used (al- 
though the changes of u affect the shape of oscillations). It proofs therefore that 
the increase of amplitudes of oscillations observed between e + e~ and pA reactions is 
caused mainly by different experimental cut-off procedures (quantified here by dif- 
ferent values of variable u used in the respective processes) applied to the measured 
P(n). In e + e~ processes, with smallest amplitudes of oscillations, the P(n) were 
measured most accurately, up to the very high multiplicities (i.e., to large values of 
the ratio z = n/(n)). The opposite situation is encountered in pA processes. This 
is the main result of our note. 



This kind of universality (even if only approximate) makes the sizes of amplitudes 
of oscillations not particularly sensitive to the dynamical details of P(n) of interest. 
Not much is left in this observable when different experiments, but with the same 
values of variable u, are compared with each other. On the other hand, the character 
of oscillations, as visualised by their frequency in the rank q of moments, remains 
in a visible way different for different types of reactions and can therefore be used 
for dynamical discrimanation between different models. For example, in Fig. 2 we 
show oscillations of H q moments obtained for the same value of u = 7 using P(n) 
in the form of MNBD as given in || |J : 

P(0) = ^ + ^" 

(1 +r 2 ) 
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p , v 1 fri\ A „ r(fc + n + j) f r 2 - ri y r 2 " 

P{n) ~ n\ \7J £ ^ r(fc + j) l^^J (4) 

where AT, fc, 7*1,2 are parameters. Referring to 0, |J for details we shall say only 
that if k = 0, the summation in eq. (fj) runs from j = 1 up to j = N and the 
resultant distribution is called MNBD. In this case parameters r x 2 are given by 
the average multiplicity (n) and second moment Ci of the corresponding Pin) as 
r i,2 = \ (p2 — 1 — -nx =F ^) (n) . For iV = 0, parameter n disapears from (||) and it 
reduces to the NBD with r% — ^u~- The parameters n and r 2 of the NMBD reflect 
now the structure of oscillations rather than their amplitudes. As one can see they 
change systematically from parameters describing e + e~ annihilation (left-top panel) 
to those typical for hadronic pp collisions (right-bottom panel) || [| • 

To summarize, we stress again that the magnitude of observed oscillations of H q 
moments of multiplicity distributions P(n) reflect essentially our ability to measure, 
in a given reaction, large multiplicities. When analysing data using the same value 
of our universal cut-off parameter u one gets comparable values of amplitudes for all 
reactions of interest. It means that this quantity is not sensitive to dynamical details 
of reaction. The shape of oscillations remains, however, sensitive to such details. It 
can therefore be used to extract a new dynamical information from different mul- 
tiplicity distributions (when compared at the same values of the cut-off parameter u) . 

The separate issue is the problem of oscillations in heavy ion collisions AB, which 
we should now briefly address for the sake of the completeness of presentation. They 
do not share the property discussed above. The reason is the following. As was 
already mentioned, in the collisions of two nuclei, A and B, the nuclear geometry 
is the main factor responsible for the shape and properties of the corresponding 
multiparticle distribution of produced secondaries P(n) [IT]. This fact is crucial 



in generating oscillations in the respective cumulant moments. To show it on some 
example let us first write the typical corresponding multiplicity distribution for A+B 
collision: 

(Hot fJ- / M \ 

P(n) = £ p{n) J] Piirn) dln-^nA. (5) 

Ai=l i=l V i=l } 

It contains two ingredients: distribution p(fi) of the number of emitting sources 
/i and the respective "elementary" multiplicity distributions of particles produced 
from such sources, Pj(rij). The emitting sources can be, for example, understood 
as in [jn|. Their distribution can be calculated in the same way as in [10[ [11, [TJ 



In the example below we have used a simple Monte Carlo code in which two col- 
liding nuclei consisting of A and B nucleons, respectively, collide with each other. 
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Nucleons are distributed in nucleus according to a standard Saxon- Woods (SW) 
distribution (with diffusiveness 0.49 fm for S and 0.545 fm for Pb nuclei and corre- 
sponding nuclear radii given by the formula: r[fm] = 1.12A 1 / 3 - 0.86/A 1 / 3 ). They 
collide with each other with probability given by their (total inelastic) cross section 
a = 32 mb. This provides us with p(fi). On the other hand Pj(rij) has been taken 
again from the MNBD fits to elementary collisions performed in |3], |3J]). In nuclear 
collisions two distinct classes of events occur and must be treated separately: central 
and minimum bias collisions. In our case central collisions were chosen as 1% of the 
collisions with smallest impact parameter. In Fig. 3 we show results for moments 
H q of P(n) from eq. @ calculated for S + S (left panels) and Pb + Pb (right pan- 
els) minimum bias (upper panels) and central (lower panels) collision for the same 
values of the variable u = 5. Notice that the magnitude of amplitude of oscillations 
(especially for central collisions) remains different from that in the corresponding 
panels of Fig. 1. It means that in this case there is no such universality as in the 
previously discussed reactions. On the other hand, however, minimum bias colli- 
sions are distinctively different from the central ones, which show only very small 
oscillations. The patterns shown apparently depend only weakly on the choice of 
the colliding nuclei (i.e., on the parameters of the Monte Carlo producing p(fi)). 



To understand better results presented in Fig. 3 one should realise that central 
collisions result in a large number of elementary collisions, i.e., in a large number of 
emitting sources fi in each event. Therefore, because of central limit theorem, irre- 
spectively of details of elementary collisions P(n) must have a gaussian-like shape. 
We can parametrize it as: P{n) = P ■ exp ^— '■""""^ On the other hand, the min- 
imum bias collisions result in p(fi) of the box-like, or Saxon- Woods (SW)-like shape 
and such will be also resultant P(n): P(n) = 1 + e (f-» o) / a - m Fig- 4 we show, as 
illustration, some typical examples of oscillation patterns emerging from both types 
of distributions. Notice that whereas we essentially observe no oscillations in the 
case of gaussian P{n) (or, if at all, they do start at large q), we see strong oscillations 
for the SW P(n). They are caused in this case by the box-like shape of the SW 
distribution, which is best demonstrated by the fact that they gradually vanish with 
the increasing difussiveness of SW used, i.e., with the increasing values of parameter 



a ]nj. It should be pointed here that results presented in Fig. 4 were obtained 
without additional truncation in multiplicity, i.e., in (|3[) n tot = oo. All oscillations 
present there are thus entirely of different origin than the simple truncation of P(n). 
They are governed by the geometrical parameter a and by the level of observability 
of the total P(n). 
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Summarizing, we have demonstrated (approximate) universality of amplitudes of 
oscillations of cumulant moments when compared at the same values of the variable 
u as defined in (Q). It shows up for a range of reactions from e + e~ annihilation 
processes, via hadronic pp collisions, to pA reactions. The latter start to show influ- 
ence of the geometry of collision process, which entirely dominates the truly nuclear 
collisions. 
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Figure captions 



Fig. 1 Examples of H q moments of multiplicity distributions for the e + e~ annihila- 
tion (upper panels), pp reactions (middle pannels) and pA reactions (bottom 
panels) for two chosen values of the parameter u: u = 5 (left panels) and u = 8 
(right pannels). (The P(n) data are the same as in Jj], 0, [|, f|, |10j). 

Fig. 2 The H q moments obtained from the MNBD for P(n) for different values of its 
characteristic parameters r\ and r 2 (cf. |3|, |4]]). The upper- left panel corre- 
sponds to e + e~ and bottom-right one to pp reactions, respectively. The value 
of parameter u = 7 reamains all time the same. Notice the gradual change 
of frequency of oscillations whereas their amplitudes remain essentially of the 
same order of magnitude. 

Fig. 3 The H q moments calculated for S + S (left panels) and Pb + Pb (right panels) 
collisions of the minimum bias (upper panels) and central (lower panels) type. 
In both cases u = 5. 

Fig. 4 Examples of oscillation patterns for gaussian-like (upper panels) and SW-like 
(lower panels) shapes of multiplicity distributions P(n), see text for details. 
Short dash, long dash and full lines correspond to parameter a equal to a = 
2, 20 , 80 for gaussian-like distributions and to a = 99, 50, 10, 0.001 for 
SW-like distributions; in both cases no = 400. 
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Figure 1 




Figure 2 




Fig. 3 
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Fig. 4 
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